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In this work we study quantum phase transitions and persistent currents in capacitively coupled 
one-dimensional Josephson-junction arrays. We will focus particularly on the roles of exciton-like 
pairs in the strong coupling limit in the presence of external gate charges and magnetic fluxes. We 
use the numerical density-matrix renormalization group method for the study in the full range of 
values of gate charge and magnetic flux. To clarify the various effects, we report the pair correlation 
functions and the exciton densities as welU as the persistent current. 
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I. INTRODUCTION 

Systems of Josephson junctions between small super- 
conducting grains have been attracting considerable in- 
terest for more than two decades. One of the main attrac- 
tive features is that they exhibit manifestations of various 
phenomena in diverse fields of condensed-matter physics. 
A popular example in contemporary mesoscopic physics 
is the Coulomb blockade effect and single charge (electron 
or Cooper pair) tunneling.^'i Persistent current, another 
hot topic in mesoscopic physics, can also be embodied 
in Josephson-junction systems. Since the superconduct- 
ing coherence is easily maintained over a macroscopic 
length scale, a "necklace" of Josephson junctions (i.e., 
one-dimensional periodic array of Josephson junctions) 
may be a good testbed for persistent currentsi^ More- 
over, charge fluctuations present in such systems may in- 
duce quantum phase transitions^'^'SiS providing a pro- 
totype model for the noble many-body phenomena in 
strongly correlated electron systemsi^ Another important 
and appealing feature of the systems is the experimental 
tunability: They not only make mesoscopic devices on 
their own^-^" but also allow us to test and understand 
otherwise very subtle points of interacting many-particle 
systems^ which is important from a fundamental point 
of view. 

Here we consider a particular geometry of Josephson- 
junction systems: a ladder of two capacitively coupled 
one-dimensional (ID) Josephson-junction arrays. In the 
Coulomb blockade regime, a single charge cannot tunnel 
across the junction since it is energetically unfavorable. 
Transport is therefore dominated by more complex ele- 
mentary processes that involve several charge-tunneling 
events. For the particular type of coupling through large 
inter-array capacitances, the relevant elementary process 
in the absence of the gate charge consists of cotunnel- 



ing of bound pairs of excess and deficit charges, which 
we call "particles" (excess charges) and "holes" (deficit 
charges), respectively It was first demonstrated on ca- 
pacitively coupled normal-metal tunneling junctionsi^ii^ 
and later on superconducting junction arrays^^^ In the 
presence of gate voltage applied between the electrode 
islands and the substrate, the particle-hole symmetry is 
broken and the particle-hole pair no longer makes the 
lowest charging-energy configuration. For example, when 
the particle-hole symmetry is broken maximally (corre- 
sponding to the gate charge given by one half of the el- 
ementary charge 2e), the transport is governed by co- 
tunneling of particle-void pairs (with the void denoting 
the absence of any excess or deficit Cooper pair)ii^ It is 
noted that these particle-hole pairs or particle- void pairs 
are reminiscent of excitons, i.e., bound states of a band 
electron and a hole, in solids. In the previous work,— 
quantum phase transitions induced by the cotunneling 
of particle-hole pairs and particle-void pairs near the 
particle-hole symmetry line and the maximal-frustration 
line, respectively, have been studied by means of pertur- 
bative methods. However, properties of the transport or 
phase transitions in between have not been studied. 

Effects of an external magnetic flux threading the 
loop of a ladder of two capacitively coupled Josephson- 
junction necklaces (CCJJNs) (see Fig. ^ are even more 
sophisticated since the objects involved in the persistent 
current are not single charges. Unlike most studies of 
the persistent current (or equivalently, the underlying 
Aharonov-Bohm effect), which focus on single-charged 
particles, recent researches into a nano-structure with 
non-simply-connected geometry have demonstrated^^ 
that excitons can contribute to persistent currents, in 
spite of their charge neutrality. The nonvanishing per- 
sistent currents in the system, is attributed to the fi- 
nite probability of breaking and recombination of an ex- 
citon via intermediate single particle/hole states. It is 
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thus quite intriguing to investigate persistent currents in 
CCJJNs, where cotunnehng of the particles and holes of 
Cooper pairs dominates the transport phenomena. Ad- 
ditional advantage of the CCJJNs is that the particle- 
hole or particle- void pairs are stable while the excitons in 
semiconductor nano-rings usually have a finite and short 
life time. Notice further that the CCJJNs are already 
within the reach of experimental realization. IS, 

In this work we study quantum phase transitions and 
persistent currents in a ladder of CCJJNs. We focus 
particularly on the roles of "excitons" in the presence 
of the charge frustration due to an external gate volt- 
age and the magnetic frustration due to an external 
magnetic flux threading the necklaces. We use the nu- 
merical density-matrix renormalization group (DMRG) 
method^^ to probe the full ranges of the gate charge 
and the magnetic flux. Although we are mainly inter- 
ested in the strong-coupling limit, we will consider for 
comparison both the two limiting cases: decoupled and 
strongly coupled cases. In the limit of strong coupling, we 
identify two different superfluid phases, characterized by 
condensation of either particle-hole pairs or particle- void 
pairs, depending on the gate charge. In order to dis- 
close properties of the superfluid phases and formation 
of excitons explicitly, we measure the pair correlation 
function and the exciton density. The behavior of the 
persistent current calculated for small systems reveal the 
transport via the separation and recombination process 
for small Josephson energies. At larger Josephson ener- 
gies, however, the transport is governed by mixing of low- 
lying charge states with higher-energy states. Finally, we 
propose an experimentally realizable system to demon- 
strate the cotunneling process of excitons. The interme- 
diate coupling regime is also interesting and more feasible 
experimentallyii^ Unfortunately, however, the numerical 
DMRG study in this case is beyond the current comput- 
ing power, requiring far more memory than available. We 
thus leave the intermediate region for the future study. 

The remaining part of this paper is organized as fol- 
lows: In Sec.|n]wc first describe the model Hamiltonians 
and discuss qualitatively the relevant low-energy charge 
states. Quantum phase transitions in a single Josephson- 
junction necklace and capacitively coupled Josephson- 
junction necklaces are examined in Sec. lIIII Section llVl is 
devoted to the investigation of the persistent currents in 
the system, revealing the AB effect of excitons. Finally, 
we summarize the main results in Sec. 



II. MODEL 

We consider two ID periodic arrays, which we call 
necklaces, of N superconducting grains as shown in 
Fig. n Any two nearest-neighboring grains on one neck- 
lace form a Josephson junction of coupling strength Ej. 
The two necklaces are coupled with each other via capac- 
itance Cj between corresponding grains, to form a "lad- 
der" . Uniform gate voltage Vg is applied to each grain 




(b) 

FIG. 1: Schematic diagrams of (a) a single Josephson-junction 
necklace (from above) and (b) two capacitively coupled neck- 
laces (from the diagonal direction). In (b) each thick ring 
represents the single necklace depicted in (a). 

through its self-capacitance Co, inducing gate charge 
Q = CoVg on each grain. For convenience, we measure 
the charge in units of 2e, and write Q = 2eng. In ad- 
dition, a transverse magnetic flux threads each necklace. 
Such a system is described by the Hamiltonian 

-EjJ2cosi^i-^'.+i-A.), (1) 

Ix 

where the number of the Cooper pairs and the phase 
(j)^. of the superconducting order parameter at site x on 
the ^th necklace {1 — 1,2) are quantum- mechanically con- 
jugate variables: (/'J^/] = iSi^i'Sx.x'- The bond angle 
Ax is given by the line integral of the vector potential A 
introduced by the applied magnetic field: 

where / denotes the total flux in units of the flux quan- 
tum $0 = 2TThc/2e. Assuming that junction capacitances 
are negligible, we write the capacitance matrix C'ix,i'x' hi 
the form.^^ 

Cix,i'x' — [CqSi^i' + Ci{2Sij' — 1)] Sx.x' = Cii'Sx.x' (3) 
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FIG. 2: Energy levels of the charging energy part in Eq. Q 
and corresponding charge configurations near the particle- 
hole symmetry line. Filled and empty circles denote particles 
and holes, respectively; paired (upper and lower) solid lines 
represent the two coupled arrays, the couplings between which 
are illustrated by the dashed lines. The low-lying energy lev- 
els satisfying = are well separated by a large amount of 
energy (of the order of Eo) from those with ^ 0. 



FIG. 3: Energy levels and corresponding charge configura- 
tions near the maximal-frustration line. Note that the ground 
state is two-fold degenerate per site. 



strong coupling regime. To examine the charge configu- 
rations in the low-lying states, it is convenient to rewrite 
the charging energy part in the Hamiltonian Q): 



and also define charging energy scales Eq = e^/2Co and 
Ej = 6^/2(7/, associated with the corresponding capac- 
itances. Notice that when rig = the Hamiltonian in 
Eq. is symmetric with respect to the particle-like (ex- 
cess Cooper pairs) excitations and hole-like (deficit pairs) 
ones. On the other hand, charges on each grain are max- 
imally frustrated when Ug = 1/2. For later use, we thus 
name the lines corresponding to rig = and Ug — 1/2 the 
particle-hole symmetry line and the maximal-frustration 
line, respectively, in the phase diagram. 

For the DMRG analysis, we represent the Hamiltonian 
in Eq. Q in the boson number basis. Based on the com- 
mutation relation between the number and the phase 
(/)^, we identify the boson creation operator b'-J at site x 
on necklace I with e*"^^ . In terms of the boson operators, 
we thus obtain the Bose-Hubbard Hamiltonian 

ll'x Ix 

-§E(«"'"'^''^^'^Ui + H> (4) 

Ix 

where n is the average boson number per site. Note that 
in the quantum phase model the Josephson energy term 
is independent of the number fluctuations, while the cor- 
responding (hopping) term is not in the Bose-Hubbard 
model. To alleviate the effects of number fluctuations in 
the Bose-Hubbard model, we consider the case that the 
average boson number n per site is largeiS: Throughout 
this study we set n to be 10000. 

Capacitive coupling between necklaces drastically af- 
fects the low-lying charge excitations, especially, in the 



J^BH = Uo E« - 2%)' + ^0 E«)' 

X X 

Ix 

where we have defined new energy scales Uq = 2Eq and 
Vo = 2i?o/(l + 2C//Co) and charge variables = 
— ri) ± {n2^x — n). Note that n+ represents the 
total number of (excess) Cooper pairs on the xth rung 
of the ladder. In the regime of concern (C/ ^ Cq, i.e., 
Ej <C Eq), we have Uq » Vq and Uq :§> Ej, and the 
charge configurations satisfying — 2ng ~ are thus 
strongly favored. In such a charge configuration /2 
corresponds to the number of excitons (particle-hole or 
particle- void pairs, see below). 

The representation in Eq. ISJ of the Hamiltonian allows 
us to distinguish clearly the two interesting regions from 
each other: near the particle-hole symmetry line {ug = 0) 
and near the maximal-frustration line (rig = 1/2), as 
one can observe from the energy spectra of the charg- 
ing energy part illustrated in Figs. El and |3| for the two 
regimes, respectively. Near the maximal-frustration line, 
the charge configurations that do not satisfy the condi- 
tion ri+ = 1 (for all x) have a huge excitation gap of 
the order of Eq. Furthermore, the ground states of the 
charging energy part, separated from the excited states 
by a gap of the order of Ej, have two-fold degeneracy 
for each x, corresponding to n~ = ±1. Near the particle- 
hole symmetry line, on the other hand, low-energy charge 
configurations should satisfy the condition n+ = for all 
X. Unlike the former case, the ground state of the charg- 
ing energy is non-degenerate and forms a Mott insulator 
characterized by rii J. = n2^x = for all x. As Ej is 
turned on, the ground state is mixed with the states with 
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n~ = ±2. In the intermediate region (0 < Ug < 1/2), 
these two kinds of energy spectra are interleaved to form 
a complex shape of the energy levels. 

III. QUANTUM PHASE TRANSITIONS 

The competition between charge order and phase co- 
herence gives rise to quantum fluctuations and quantum 
phase transitions at zero temperature. For large charg- 
ing energy {Eq » Ej), the bosons become localized and 
the system is in the Mott insulator phase with integer 
density. On the contrary, for large hopping energy or 
Josephson energy {Ej 3> i?o), coherence of the phases (j)x 
dominates over the system and the superfluid (SF) region 
with delocahzed bosons is observed. The properties and 
universality classes of the phase transitions, however, de- 
pend strongly on the coupling strength C//Co as well as 
the chemical potential fi = SEoing + h). The charge frus- 
tration Hg may be restricted to the range [1, 1/2] since the 
Hamiltonian in Eq. ^ is periodic in Ug with period unity 
and has reflection symmetry about the ng = 1 (or any 
integer) line. In the foUowings we investigate two limit- 
ing cases: the decoupled case (C/ = 0) and the strongly 
coupled one (C/ 3> Co). 

A. Single Josephson-junction necklace 

The phase transition in a single Josephson-junction ar- 
ray has been studied quite extensively and it has been 
found that its nature depends crucially on the gate volt- 
age. In the presence of nonzero gate voltage {ug ^ 0), 
the density of the system changes as the phase bound- 
ary is crossed from the incompressible insulator to the 
compressible superfluid. The transition can thus be lo- 
cated at the point where in the thermodynamic limit the 
density of the ground state becomes different from one 
of the insulator ground state as the Josephson energy is 
increased. On the other hand, in the particle-hole sym- 
metry line {ug — 0), the density remains to be an integer 
at the phase transition. Therefore, in this case the phase 
boundary is determined by the single-particle excitation 
gap. This is possible because in the superfluid phase the 
ground state is a superposition of states with different bo- 
son numbers, the energy gap between the ground state 
and the states with additional particles, which is finite in 
the insulating phase, vanishes in the superfluid phase. 

Since the Hamiltonian conserves total charge number, 
the DMRG algorithm can be set up to target states with 
given total excess number M of bosons. We thus obtain 
the phase diagram of the system by comparing energies of 
the ground states with different boson numbers: the en- 
ergy Em=o of the insulator ground state with zero excess 
boson density ((n^) = fi) and the energy Em=i of the 
eigenstate with an additional particle upon the ground 
state. Through the linear extrapolation of the energy gap 
Em=i —Em=q for finite system size N = 64, 128, and 256, 
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FIG. 4: Phase diagram for a single ID Josephson array. The 
phase boundary separates the Mott insulator (Ml) with zero 
excess Cooper pair density from the superfluid phase (SF). 
The Josephson energy Ej is expressed in units of 8Eo and 
the error bars in the Ej direction (for given Ug) are smaller 
than the square symbols. The line is merely a guide to the 
eye. 



we have estimated the gap in the thermodynamic limit 
and located the phase boundary at the point where the 
gap is zero. For high numerical accuracy and access to 
large systems, the finite-size DMRG algorithm and open 
boundary conditions have been used. During the DMRG 
process, the boson number at each site is truncated to 
be less than six and the discarded weight is set to be 
less than 10~^, giving rise to negligible errors in the gap 
energyiSHi The magnetic frustration / is set equal to zero 
because it can be gauged away and becomes irrelevant in 
the thermodynamic limit. 

Figure ^ displays the resulting phase diagram for the 
quantum phase model on the Ej-Ug plane, in the range 
< JT-g < 1/2 and < Ej < SEq. For convenience, here 
and in all subsequent figures, the energy is expressed in 
units of SEq. The phase diagram, where the Mott insu- 
lator region with zero excess boson density is separated 
from the compressible superfluid, is in good agreement 
with those obtained via the perturbative expansionSi and 
the quantum Monte Carlo method.-— On the particle- 
hole symmetry line the quantum phase model is mapped 
exactly to the (l+l)-dimcnsional XY model, predicting 
a Beresenskii-Kousterlitz-Thouless (BKT) transitions^ 
driven purely by phase fluctuations. The sharp-pointed 
shape of the insulating region near the symmetry line re- 
flects the slowness in closing the energy gap in the BKT 
transitioniSS In case that the particle-hole symmetry is 
broken (away from the symmetry line) , no such slowness 
is found and the commensurate-incommensurate transi- 
tion belongs to a universality class other than that of the 
XY model, with different critical exponents^^ and RG 
characteristics!^ 



B. Strongly coupled Josephson-junction necklaces 

In the strong coupling limit, the low-energy charg- 
ing states relevant to the phase transition are the the 



5 



particle-hole pairs (with n+ = and n~ = ±2) and the 
particle- void pairs (with = 1 and — il)- For 
small hopping strength, these excitons are localized and 
the system is in the Mott insulator phase. As Ej in- 
creases, the phase boundary is crossed from the insulator 
to the superfluid which, in this case, originates from con- 
densation of the excitons. Accordingly, as in the case of a 
single Josephson-j unction necklace, the transition can be 
located as one tracks the energy taken to add an exciton 
to the insulator: At the phase boundary this energy van- 
ishes in the thermodynamic limit. Which kind of exciton 
between the particle-hole pair and the particle-void pair 
is relevant depends on the charge frustration Ug. Near the 
particle-hole symmetry line {ug « 0), the particle-hole 
pairs are energetically favorable and governs the phase 
transition. As Ug is increased and in the presence of the 
Josephson tunneling, in contrast, the particle-void pairs 
begin to be dominant faster than the particlc-hole pairs, 
which will be shown below. 

In the DMRG procedure we have associated the target 
state with a pair of total excess boson numbers (Mi, M2) 
on the two arrays by utilizing the boson number conser- 
vation. In order to locate the phase boundary, we have 
calculated the energy E(^Mi,M2) of three kinds of eigen- 
state: the insulator ground state with (Mi,M2) = (0,0) 
and the states with additional particle-hole and particle- 
void pairs upon the ground state, labeled by (Mi, M2) = 
(1,-1) and (1,0), respectively. We have extrapolated 
the energy gaps £'(i,--i) - ^^(co) and £'(i_o) - £'(0,0) for 
finite system size N =16, 32, and 64 to locate the transi- 
tion points where the gaps vanish in the thermodynamic 
limit. As in the case of a single Josephson necklace, we 
have employed the finite-size DMRG algorithm, imposing 
open boundary conditions. 

Figure [3 shows the energy gaps as functions of the 
Josephson energy at various charge frustrations, in the 
system with C7/C0 = 100 and N = 32. From Fig. E^a) 
we observe that the excitation energy i?(i _i) — £^(o,o) 
for different charge frustrations collapses into one curve, 
which also happens at other system sizes. This indicates 
that the critical Josephson energy at the transition driven 
by the particle- hole pairs does not depend on ng. On 
the other hand, the energy gap -E(i,o) ~ -£'(0.0) decreases 
almost linearly with the increase of Ug and Ej, as shown 
in Fig. Elb). The larger ng is, the smaller the Josephson 
energy Ej at which the energy gap vanishes becomes. 
For Ug > 0.14, the critical value of Ej become even less 
than that for the particle- hole pairs. 

The resulting phase diagrams for strongly coupled ar- 
rays with Ci/Cq = 100 and 200 are exhibited in Fig. El 
(a) and (b), respectively. Based on the dominant trans- 
port mechanism, one can distinguish three regions in the 
superfiuid phase: ESi, ES2, and SFUB. In region ESi the 
transport is driven mainly by the excitons of particle-hole 
pairs; in ES2 it is driven by particle-void pairs. In region 
SFUB, on the other hand, single-particle processes dom- 
inate the transport in the system. Such superfiuid of un- 
paired bosons (SFUB) is to be observed at Ej/SEq ^ 1, 
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FIG. 5: Energy gaps (a) -Bfo.o) and (b) £"(1,0) -^(o.o) . 

Both are taken for Ci/Co — 100 at the system size A'' = 
32. Each symbol corresponds to a different value of charge 
frustration: (a) Ug = 0(0), O.l(o), 0.2(A), and 0.24(0); (b) 
Ug = 0.14(0), 0.16(B), 0.18(o), 0.2(.), 0.22(A), and 0.24(a). 
The energies on both axes are expressed in units of SEq. 



far to the right from regions ESi and ES2, and not shown 
in the phase diagram given by Fig. El We note that differ- 
ent transport mechanisms take over dominant roles grad- 
ually as the control parameters are changed. Therefore, 
regions ESi and ES2 in Fig. El should not correspond to 
truly distinct phases. 

Previous studies on mapping of the system at the 
particle- hold symmetry line to (H-l)-dimensional sys- 
tem of classical vortices2SiS& insisted that the system is 
effectively described by a two-dimensional XY model 
and exhibits a BKT transition at the critical Joseph- 
son energy Ej/8Eo\c = 4i^|^T(l + + 2Ci/Co)-^ « 
2Kgp.rp{Co/Ci), where Kbkt ~ 0.748 is the critical cou- 
pling strength for the standard XY model. Our data, 
though being unable to discern nature of the transition, 
show that the critical Josephson energy is inversely pro- 
portional to yJCijC^ instead of C/ /Co, apparently favor- 
ing against the BKT transition. This result is quite rea- 
sonable in view of the fact that the cotunncling process 
of particle-hole pairs via an intermediate virtual state 
happens with the probability proportional to Ej/EqEj, 



leading to Ej/Eq\c oc ^JCijC^. In addition, the nonzero 
charge frustration does not change the properties of the 
phase transition abruptly, in contrast to the case of a 
single array. Instead, the transition point as well as the 
qualitative properties is preserved up to Ug « 0.135 for 
C//Co = 100 and to Ug w 0.165 for C//Co = 200; there 
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FIG. 6: Phase diagrams for strongly coupled ID Josephson 
arrays for (a) Ci/Co = 100 and (b) Ci/Co = 200. Displayed 
are the phase boundaries between the Mott insulator (MI) 
phase and the superfluid phase (ESi and ES2). Regions ESi 
and ES2 in the superfluid phase are distinguished by the dom- 
inant transport mechanism (see the text). The transitions 
across the boundaries indicated by filled and empty symbols 
are driven by the particle-hole pairs and by the particle- void 
pairs, respectively. The lines are merely guides to eyes. 
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FIG. 7: Pair correlation Cpair versus the Josephson en- 
ergy, (a) For Ug = 0, the pair correlation functions in 
the insulator ground state (empty symbols) and in the state 
with one particle-hole pair added upon the ground state 
(filled symbols) are plotted, (b) shows the correlations in 
the ground state in the presence of the gate voltage Ug = 
0.2(n), 0.3(0), 0.4(A), and 0.5(v). Here we set Ci/Co = 100 
and the system size N = 8, and lines are guides to eyes. 



is no increase in the critical value of Ej as predicted in 
Ref. Q. Since our model neglects the junction capaci- 
tance on each necklace, for 

each site has two-fold degenerate ground states (n+ — 1 
and n~ = ±1) of the charging energy. Accordingly, the 
Josephson energy of any strength brings about charge 
fluctuations to drive the system into the superfluid phase. 
Indeed Figure shows that the MI phase ceases to exist 
for Ug > n*, regardless of Ej. With nonzero junction 
capacitance, the degeneracy is expected to be broken, 
generating another insulating phase: the charge-density 
wave (CDW) phase. The perturbative study^'* has found 
that as the Josephson energy is increased the system goes 
from the CDW insulator to the Luttinger liquid phase. 

To witness the activity of the excitons in the phase 
transition, we have measured the pair correlation func- 
tion defined to be 

X 

together with the exciton density Pq of the particle-hole 



pairs and P±i of the particle-void pairs: 

■fo = -^^{Snl^+nlfi - Sni..oSnl,o) (8) 

X 

^ ^E('^"i+"^^«+) (n+^O)- (9) 

X 

The pair correlation function assumes zero if there is no 
correlation between the boson numbers on the two arrays; 
a particle-hole or particle-void pair at every site on the 
average contributes to Cpair by —1 or —1/4. 

FigureHshows that the pair correlation is negative and 
monotonically decreases with Ej, which indicates that 
larger hopping strength makes more excitons come into 
the system. For Ug — [see Fig. Ufa)], the two pair cor- 
relation functions, one for the insulator ground state and 
the other for the state with an additional particle-hole 
pair, approach each other and collapse at Ej/SEq > 0.07, 
giving another evidence for the condensation of excitons. 
For Ug = 0.2, the correlation changes abruptly at the 
phase transition, as shown in Fig.[7Jb). For Ug > 0.3, we 
have Cpair « — 1/4 for small values of Ej, which reflects 
the contribution of particle- void pairs. Although the cor- 
relations for Ug > 0.3 become similar at large values of 
Ej, they do not coincide at intermediate values, implying 
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FIG. 8: Exciton densities (a) Pq and (b) Pi versus the Joseph- 
son energy for gate voltage Ug — 0(0), 0.2(B), 0.3(o), 0.4(»), 
and 0.5(A). Other parameters are the same as those in Fig.Q 
The pairs of numbers (Mi, M2) represent the total excess bo- 
son numbers of the ground state at given parameters. 



FIG. 9: Persistent current / and ground state energy E as 
functions of the flux / in (a) the insulating phase {Ej/8Eq = 
0.64) and (b) the superfiuid phase {Ej/8Eo — 1) along the 
particle- hole symmetry line {ug — 0). 



IV. PERSISTENT CURRENT 



the difference in the history of pair generation. 

The behavior of the exciton density with Ej, displayed 
in Fig. |S1 manifests more clearly the formation of the 
excitons. The particle-hole pair density Pq grows with 
Ej [see Fig.lHIa)], except for the case of Ug ~ and large 
Ej. On the other hand, Fig.|H{b) shows that the density 
Pi of the particle-void pairs with nl. + ri^ = 1 has two 
kinds of tendencies: For Ug < n* the density Pi increases 
with Ej; otherwise it decreases. The steep changes in 
both densities Pq and Pi happens when the total excess 
boson number M — Mi + M2 of the ground state is 
altered. Such behaviors of Pq and Pi reveal that more 
than one kinds of excitons proliferate in the system as 
Ej is increased beyond its critical value. With large Ej, 
the kinetic energy gain due to the Josephson tunneling 
term compensates for the charging energy gap between 
different kinds of excitons. 

We close this section with a comment about the pair 
correlation in the limit Ej/8Eo ^ 1, which is beyond 
our current computational power. Our data shows no 
indication for the decrease of the pair correlation with 
Ej raised. However, when the Josephson energy is large 
enough for the single processes of unpaired particles to 
prevail, the pair correlation may eventually approach zero 
again. 



In this section we consider the persistent current along 
the necklaces, induced by the threading external mag- 
netic field. Since tunneling of Cooper pairs between neck- 
laces is not permitted, the persistent current carried by 
each necklace is given by the derivative of the energy with 
respect to the magnetic flux 

(10) 

which is simply the supercurrent through the Joseph- 
son junctions. The current in the system is thus given 
by the imaginary part of (exp — ^^+1 — A^)] ), the 
real part of which describes the gauge-invariant phase 
correlation function between nearest neighboring grains. 
Since the current is periodic in / with period unity and 
an odd function of the flux /, it is sufficient to calculate 
the current in the range / e [0, 1/2]. As in the previous 
section, we focus on two extreme cases: the decoupled 
case (C/ — 0) and strongly coupled one {Cj/Cn ;» 1). 

A. Single Josephson-junction necklace 

We have calculated the persistent current in a flnite- 
size system {N = 40) under periodic boundary condi- 
tions. The persistent current is evaluated in the ground 
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FIG. 10: Persistent current / versus the system size A'^ in 
(a) the insulating phase [Ej = 0.25(B), 0.36(»), and 0.49(a) 

in units of 8Eo] and (b) the superfiuid phase [Ej — 0.81(n), 
l(o), and 1.44(a) in units of 8Eo] along the particle-hole sym- 
metry line {ug = 0). Solid and dotted lines represent the 
exponential function ae~^^ and the algebraic function c/N, 
respectively, where the constants a, b, and c are obtained via 
a fitting algorithm. 



state, namely, the lowest-energy state, which is found by 
varying the total excess boson number M at given Ug and 
Ej, and expressed in units of eEj/Nh in all subsequent 
figures. Figures E^a) and l^b) show the dependence of 
the persistent current (in units of eEj/Nh) and of the 
ground state energy (in units of 8Eq) on the flux / in 
both the insulating phase and in the superfiuid phase, 
respectively, without the gate charge (rig = 0). For small 
Ej, the current depends sinusoidally on /, whereas it 
has a saw-tooth shape in the superfiuid phase. Such be- 
havior of the persistent current is well known in the two 
extreme one-dimensional electron models: In the tight- 
binding model with the lattice potential energy domi- 
nant over the kinetic energy, the single-particle energy is 
given by a cosine function of the flux /, giving rise to 
sinusoidal dependence of the current on /. On the other 
hand, the free electron model on a ring, where the energy 
is quadratic in /, has the persistent current linear in / 
and of the saw-tooth shape. In our model Cooper pairs 
take the role of the electrons and in analogy we infer that 
the saw-tooth dependence in the persistent current indi- 
cates the emergence of the superconductivity over the 
system, where the Cooper pairs can freely move around. 
Our data for the dependence of the current on the system 
size, shown in Fig. ^1 also leads to the same interpreta- 



0.03 





4 


1/ 

1 'i 




A-' 


• ■ 












□ 






n=Q.l 


■ 


• 




n=Q.2 






■ 1/ 


n=0.3 - 




: • 




n=OA 


- - -A- - - 










A ■ O 









0.5 1 1.5 

Ej 

FIG. 11: Persistent current versus the Josephson energy for 
various gate voltages. In each case the current, scaled in units 
of eEj/Nh, makes clear the contribution of the correlations 
between nearest-neighboring grains. All the currents are cal- 
culated at / = 1/4. 



tion: In the insulating phase the bosons arc localized at 
sites so that the probability for a boson to circle around 
the ring and to return to its starting position is propor- 
tional to , where t is the hopping probability between 
nearest neighbors. This gives the current decaying expo- 
nentially with the system size N [see Fig. llOf a)]. On the 
other hand, in the superfiuid phase the wave function of 
the boson is extended and the hopping probability over 
the system does not depend on the system size. Instead, 
since the energy itself is quadratic in the system size, the 
persistent current follows a power-law with respect to the 
system size, as shown in Fig. llOf b). Hence our data for 
the persistent current are fully consistent with the phase 
transition explained in Sec. IIII Al 

We exhibit the dependence of the persistent current on 
the Josephson energy at / = 1/4 and various gate volt- 
ages in Fig. ^2 The persistent current scaled by Ej 
is negligibly small in the insulating region, then rises 
rapidly near the transition point, and increases only 
marginally in the superfiuid phase. For Ug = 0, the cur- 
rent shows finite-size effects, gradually increasing quite 
before the transition point. On the contrary, in the pres- 
ence of nonzero gate voltage, the current increases very 
sharply at the phase boundary even in the small-size sys- 
tem, which is attributed to the abrupt change in the total 
boson number or the density of the ground state at the 
transition point. Deep in the superfiuid phase, on the 
other hand, the persistent current becomes independent 
of the gate voltage. 



B. Strongly coupled Josephson-junction necklaces 



From the observation in Sec. IIII Bl it is evident that in 
the strongly coupling limit the excitons such as particle- 
hole and particle-void pairs play dominant roles in the 
transport. In the picture of the lowest-order cotunneling 
processes illustrated in Fig. 1121 however, such pairs are 
tightly bound throughout the transport process. Accord- 
ingly, the current induced in one necklace is accompanied 
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(b) 

FIG. 12: Typical cotunneling processes relevant (a) near 
the particle-hole symmetry line and (b) near the maximal- 
frustration line. Such cotunneling results in the current mir- 
ror effects. 




FIG. 13: Persistent current I\ (along one necklace) versus the 
Josephson energy Ej in the strongly coupled system of size 
iV = 8, for / = 1/4 and gate voltage Ug = 0(n), 0.2(o), and 
0.5(A). Lines are merely guides to eyes. 

by the secondary current in the other necklace, with the 
same magnitude but in the opposite direction. On the 
other hand, in response to the magnetic flux, the charges 
in an excitonic pair tend to move in opposite directions 
since their signs are opposite (with respect to the offset 
charge ng). Therefore the current mirror effect competes 
with the influence of the magnetic flux. 

Indeed, for small values of the Josephson energy 
{Ej/8Eo < 0.07), the persistent current is quite negligi- 
ble both on the particle-hole symmetry line (ug — 0) and 
on the maximal-frustration line (rig = 1/2) (see Fig. ll3ll . 
The small but still non-zero amount of persistent current 
is induced via higher-order tunneling processes. Namely, 
the charges in an excitonic pair break up, run down the 
circumferences in the opposite directions, and recombine. 
Contributions from these processes are observable only 
in a system with a small number of sites {N = 8 for the 
data in Fig. I13|) . In other words, the current mirror ef- 
fect wins the competition. It is distinguished from the 




FIG. 14: Persistent current versus the Josephson energy on 
the particle-hole symmetry line {rig = 0) in the strongly cou- 
pled system with the coupling capacitance Ci /Co ~ 100 (□) 
and 200 (o). 



behavior of the persistent current in a single Josephson- 
junction necklace, where for = the current increases 
rapidly near the transition point and any strength of the 
Josephson energy induces rather large persistent current 
for Ug = 0.5 (see Fig. Ill() . demonstrating the action of a 
different kind of charge fluctuations in the coupled sys- 
tem. 

For larger values of the Josephson energy [Ej /SEq > 
0.07), however, a considerable amount of persistent cur- 
rent flows through the system and increases with Ej. 
It can be explained by the generation of excitations with 
higher charging energies in the presence of the Josephson 
energy. As observed in Fig.|Sl with the increased Joseph- 
son energy, more of the charge states that do not satisfy 
= (near the particle-hole symmetry line) or = 1 
(near the maximal-frustration line) are now mixed with 
the lowest charging energy states. These excitations can 
carry a flnite amount of persistent current since the signs 
of the charges in a pair are not opposite now. In short, the 
magnetic frustration wins the competition with the cur- 
rent mirror effect. It is also interesting that the persistent 
currents for rig = 0(0) and for 1/2 (A) are almost the 
same. In fact, the densities of charge excitations which 
do not satisfy the lowest-charging energy condition are 
nearly the same for the two cases since the charging en- 
ergy costs for such excitations amount to the same energy 
Uq in both cases. 

Note also that unlike these two cases the persistent 
current for Hg = 0.2 (o) increases sharply at the tran- 
sition point and becomes quite larger than the one for 
Ug — or 1/2. Intermediate values of charge frustration 
(in the superfluid phase) bring about a variety of charge 
excitations in the presence of the Josephson energy and 
diminish the energy gap between the charge excitations, 
giving rise to the reduction of the current mirror effect 
and favoring independent single-charge transport. 

Figure FPil shows that the persistent current increases 
slightly as the coupling capacitance is raised. On one 
hand, a larger value of the coupling capacitance reduces 
the lowest excitation energy (~ Ej) and makes the exci- 
tons more proliferate in the system, thus increasing the 
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FIG. 15: Persistent current versus the system size on the 
particle-hole symmetry line (rig = 0) for different values of the 
Josephson energy Ej = (a) 0.045 (a), 0.08 (•), and 0.125 (a); 
(b) 0.125 (a), 0.18 (V), and 0.245 (0) (again in units of S^o)- 
The fitting curves are the same as those given in Fig. 1101 



persistent current due to the excitons. On the other 
hand, breaking of the excitons, which is crucial for induc- 
ing the persistent current, costs higher energy (~ Eq). 
These conflicting trends result in the slight increase in 
the regime of our interest {Eq ^ Ej). 

The dependence of the current on the system size also 
supports our scheme for the role of the excitons in the 
persistent current. Figure 1151 shows that, similarly to 
the case of a single necklace, in the insulating phase the 
current decays exponentially with the system size and de- 
creases inversely to the system size deep in the superfluid 
phase. Near the transition point on the side of the super- 
fluid phase (at Ej/SEq = 0.08 and 0.125), however, the 
current does display exponential dependence on the sys- 
tem size, which indicates that spatially localized objects 
participate in the generation of the current. It is another 
evidence for the virtual processes of unpaired charges or 
higher-order excitons, in the region where the low-lying 
excitons themselves are delocalized over the system. 

To reveal the cotunneling process more explicitly, we 
devise another interesting configuration that the mag- 
netic field penetrates only one necklace (/ = 1) without 
affecting the other {I = 2). Such a setup may be real- 
ized experimentally as shown in Fig.^] Notice that only 
part of the two necklaces are capacitively connected. In 
order for the persistent current to flow through uncou- 
pled grains, the Josephson coupling between those uncou- 
pled grains should be sufficiently large. In this arrange- 



FIG. 16: Schematic diagram of two coupled Josephson- 
junction necklaces. Part of the two necklaces are capacitively 
connected whereas the magnetic field threads only one of the 
two necklaces. 
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FIG. 17: Persistent currents Ji (solid lines with filled symbols) 
and I2 (dotted lines with filled symbols) when the magnetic 
field threads only the first {I — 1) necklace. The currents are 
calculated for the gate voltage Ug — (a) 0(B), 0.5(a); (b) 
0.2(»). For comparison, the persistent currents when both 
the necklaces are threaded by the magnetic field (see Fig. 
are also plotted, represented by the corresponding empty sym- 
bols. 



ment, one can observe the current mirror effect, similar 
to the case that only one chain is biased by an external 
voltageiiii^ 

Figure llTf al exhibits the current mirror effect in the 
system with the magnetic field threading only the first 
{I — 1) necklace. On the particle-hole symmetry and the 
maximal-frustration lines, the two persistent currents Ii 
and I2 along the first and the second {I — 2) necklaces, 
respectively, satisfy the relation Ii w —I2 in the range 
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Ej/8Eo < 0.07. As Ej is increased further, nonetheless, 
not only the mirror effect disappears gradually but also 
the current I2 diminishes to zero. It means that the in- 
dependent single-charge transport rather than the cotun- 
neling transport is favorable at large values of the Joseph- 
son energy. Note that the current Ii is much higher than 
the corresponding current in the system with the mag- 
netic field acting on both necklaces. Interestingly, unlike 
the previous setup, the persistent current for Ug = 1/2 is 
higher than that for Ug = 0. Whereas on the maximal- 
frustration line the system is in the superfluid state of 
the particle- void pairs even at small Ej, on the particle- 
hole symmetry line a sufficient amount of the Josephson 
energy is necessary for generating excitons, i.e., particle- 
hole pairs. For ng — 0.2, the current mirror effect is in- 
deed negligible and the increase in Ji is also very small, 
as shown in Fig. llTr b). 

V. CONCLUSION 

We have studied phase transitions and persistent cur- 
rents in a ladder of two capacitively coupled Josephson- 
junction necklaces. Emphasis has been paid on the roles 
of excitons in the presence of charge and magnetic frus- 
tration. To obtain the properties of the ground and ex- 
cited states of the system, we have utilized the DMRG 
method for arbitrary values of the gate charge and of 
the magnetic flux. Although the main interest lies in the 
strong-coupling limit between the two necklaces, we have 
studied both the uncoupled and strongly coupled cases 
for comparison. In both cases, the gate voltage brings 
about crucial effects on the properties of the system. In 
a single Josephson-junction necklace, the presence of the 
gate voltage changes rather abruptly the behavior of the 
persistent current as well as nature of the phase tran- 
sition. On the other hand, in the capacitively coupled 
Josephson-junction necklaces, such drastic change is not 



observed but the gate voltage determines the class of exci- 
tons driving the phase transition: the particle-hole pairs 
near the particle-hole symmetry line and the particle- 
void pairs near the maximal-frustration line. In the pres- 
ence of the Josephson tunneling, two different superfluid 
phases, characterized by the condensation of cither of the 
two types of excitons, have been identified, depending on 
the gate charge. The pair correlation function and the 
exciton density have provided evidence for the formation 
of such excitons. 

In the strongly coupled necklaces, the behavior of the 
persistent current is manifested by the competition be- 
tween the current mirror effect and magnetic frustration, 
which is associated with the cotunneling transport of the 
bound excitonic pairs of either particles and holes or par- 
ticles and voids. At small values of the Josephson en- 
ergy, the current mirror effect wins the competition and 
only a very small amount of persistent current is allowed 
for a finite-sized system. At large values of the Joseph- 
son energy, magnetic frustration can make use of higher- 
charging-energy states to dominate over the current mir- 
ror effect, allowing a considerable amount of persistent 
current. We have also suggested an experimentally re- 
alizable system to demonstrate the cotunneling process 
of the excitons. To out present knowledge, the only ex- 
perimental work related to the system considered here 
is that reported in Ref. ITsl Unfortunately, however, the 
large bias voltage applied to both arrays does not allow 
us to make a direct connection. In particular our DMRG 
algorithm is not suitable for such a non-equilibrium prob- 
lem. 
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